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LONG-TERM  GOALS 

The  long-term  objective  of  the  GEOCLUTTER  program  is  to  understand  the  causes  and  implications 
of  geologic  clutter  (reverberation)  in  a  geologically  well-characterized  shallow- water  environment. 

The  field  area  selected  for  the  GEOCLUTTER  program  is  the  mid-outer  continental  shelf  off  New 
Jersey,  USA.  The  New  Jersey  margin  was  chosen  for  the  GEOCLUTTER  study  because  the 
bathymetry  and  portions  of  the  shallow  subsurface  of  this  area  had  already  been  mapped  in  detail  as 
part  of  an  earlier  ONR  program  aimed  at  understanding  the  origin  of  subsurface  stratigraphy  on 
continental  margins  (STRATAFORM).  In  addition  to  multibeam  bathymetry,  ‘calibrated’  backscatter 
data  (at  95  kHz  from  the  multibeam  sonar)  was  also  collected  as  part  of  the  STRATAFORM  program. 

OBJECTIVES 

The  overall  scientific  objectives  of  the  GEOCLUTTER  program  are:  1)  to  understand,  characterize, 
and  predict  lateral  and  vertical,  naturally-occurring  heterogeneities  that  may  produce  discrete  acoustic 
returns  at  low  grazing  angles  (i.e.,  "geologic  clutter")  and  then;  2)  to  conduct  precise  acoustic 
reverberation  experiments  at  this  site  to  understand,  characterize,  and  potentially  mitigate  the  geologic 
clutter. 

APPROACH 

In  order  to  meet  these  objectives  and  to  properly  implement  acoustic  models  for  the  GEOCLUTTER 
area,  we  need  to  know,  or  predict,  the  key  acoustic  and  physical  properties  throughout  the  volume  of 
interest  (i.e.,  grain  size,  density,  sound  speed,  attenuation).  The  properties  of  the  near-surface  seafloor 
sediments  are  particularly  important.  A  possible  approach  to  this  problem  is  to  use  the  95  kHz 
multibeam  backscatter  data  collected  in  the  region,  which  may  provide  information  on  seafloor 
sediment  properties.  If  remotely-sensed  backscatter  data  can  be  used  to  infer  seafloor  sediment 
properties,  we  would  have  the  ability  to  make  quantitative  statements  about  seafloor  properties  over 
large  areas  of  the  seafloor  and  thus  the  ability  to  address  a  number  of  important  navy-related 
problems.  The  relationship  between  backscatter  and  sediment  properties  remains  ambiguous,  and  as  of 
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yet  cannot  be  used  as  a  direct  and  quantitative  predictor  of  seafloor  properties.  Attempting  to 
understand  the  relationship  between  the  multibeam  backscatter  and  the  properties  of  the  seafloor  is  the 
primary  theme  of  this  component  of  our  GEOCLUTTER  research  program. 

In  light  of  the  fact  that  we  have  not  yet  successfully  been  able  to  produce  accurate  estimates  of  seafloor 
properties  from  remotely  sensed  acoustic  data,  our  initial  proposal  fell  back  on  more  traditional  means 
of  sampling  and  laboratory  measurements  to  obtain  the  needed  seafloor  property  information  in  the 
GEOCLUTTER  area.  Given  the  coarse-grained,  sandy  nature  of  the  sediment  in  the  region  we  were 
concerned  that  laboratory  measurements  of  certain  properties  (in  particular  sound  speed  and 
attenuation)  on  core  samples  would  not  reflect  in  situ  values  as  sandy  sediments  tend  to  de-water  very 
quickly.  Thus,  the  first  phase  of  our  GEOCLUTTER  work  consisted  of  the  development  of  a  simple 
and  relatively  inexpensive  device  designed  to  measure,  in  situ,  the  spatial  variability  of  sound  speed 
and  attenuation  in  near-surface  sediments  at  the  GEOCLUTTER  site.  Our  in  situ  measurements  could 
then  be  combined  with  the  data  collected  from  cores  (by  other  investigators  -  John  Goff  from  the 
University  of  Texas  and  Chris  Sommerfield  from  the  University  of  Delaware)  as  well  as  other  acoustic 
data  (experiments  by  Charles  Holland  from  Pennsylvania  State  and  Steve  Schock  from  Florida  Atlantic 
University)  to  better  understand  the  variability  of  in  situ  sediment  physical  and  acoustic  properties  in 
the  GEOCLUTTER  area.  Subsequent  efforts  saw  the  collection  of  in  situ  physical  and  acoustic  data  at 
the  site  of  the  ONR  Martha’s  Vineyard  Mine  Burial  Experiment  (see  Mine  Burial  Annual  Report)  and 
future  work  will  involve  the  collection  of  in  situ  data  in  a  very  well  documented  region  of  Portsmouth 
Harbor,  New  Hampshire. 

WORK  COMPLETED 

Following  successful  field  work  in  support  of  the  New  Jersey  GEOCLUTTER  program  and  the  Mine 
Burial  program  off  Martha's  Vineyard,  we  recently  completed  two  small  field  experiments  in  nearby 
Portsmouth  Harbor  (PH)  and  Little  Bay  (LB),  NH.  Field  work  included  sediment  sampling  and 
analysis,  and  in-situ  measurements  with  ISSAP  ( In-situ  Sound  Speed  and  Attenuation  Probe).  In  our 
previous  field  work,  ISSAP  (see  Fig.  1)  was  configured  with  four  65  kHz  'omni-directional'  acoustic 
probes  and  'time-of-flight'  measurements  obtained  across  five  acoustic  paths.  For  this  field  work,  new 
40  kHz  'directional'  probes  were  used,  which  increased  the  frequency  range  of  the  measurements,  but 
limited  the  number  of  useable  acoustic  paths  to  two.  Prior  to  the  field  work,  the  additional  probes 
required  modifications  to  ISSAP's  transmit/receive  board  (improvements  to  the  drive  electronics  and 
increased  gain  control)  and  updates  to  the  data  acquisition  program  in  Labview. 
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Figure  1.  (a)  ISSAP  instrument  and  specifications  and  (b)  underside  view  of  ISSAP 
showing  the  configuration  of  probes. 


In  October  2003,  sediment  sampling  was  conducted  from  the  R/V  Coastal  Surveyor  and  ISSAP 
measurements  made  from  the  R/V  Gulf  Challenger.  Sediment  samples  were  obtained  with  a  Van  Veen 
grab  sampler  (University  of  Rhode  Island  (URI),  Marine  Geomechanics  Lab)  in  Portsmouth  Harbor 
(20  stations)  and  Little  Bay  (28  stations).  Sub-samples  were  taken  from  the  grab  samples  and  analyzed 
using  URI's  GeoTek  multi-sensor  core  logger  (MSCL).  Lab  measurements  including  p-wave  sound 
speed  (at  250  kHz),  gamma  ray  attenuation  (saturated  bulk  density),  magnetic  susceptibility,  and 
electrical  resistivity  were  completed  at  1  cm  intervals  on  PVC  mini-core  tubes  (5.08  cm  ID  x  8  cm 
long).  A  coarse  grain  size  analysis  (5  bins)  was  performed  on  all  grab  samples  with  an  additional 
hydrometer  analysis  completed  for  samples  with  a  mud  component  exceeding  10%.  Measurements  of 
permeability  and  s-wave  sound  speed  on  representative  samples  have  not  been  completed. 

Measurements  of  in-situ  sound  speed  and  resistivity  (porosity)  were  completed  with  ISSAP  configured 
with  two  orthogonal  matched  pairs  of  transducer  probes  operating  at  frequencies  of  40  and  65  kHz  (see 
Fig.  lb).  Measurements  were  obtained  in  both  Portsmouth  Harbor  (13  stations)  and  in  Little  Bay  (23 
stations).  In  April  2004,  ISSAP  measurements  were  repeated  in  Little  Bay  (24  stations)  using  the  same 
probe  configuration  (40  and  65  kHz  probes).  Additional  measurements  with  ISSAP  configured  using 
matched  65  and  100  kHz  probes  were  obtained  in  Little  Bay  (18  stations). 
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Figure  2.  Piscataqua  River  Watershed,  Portsmouth,  NH.  ISSAP  stations  (planned)  were  selected  to 
represent  a  range  of  sediment  types  and  are  denoted  with  a  red  symbol,  (a)  Little  Bay  field  area 
shown  (defined)  with  Simrad  EM3000  (300  kHz)  multibeam  backscatter,  (b)  Overview  showing 
relative  location  of  the  Portsmouth  Harbor  and  Little  Bay  field  areas,  and  (c)  Portsmouth  Harbor 
field  area  (defined)  shown  with  Simrad  EM3000D  (300  kHz,  dual  head)  multibeam  backscatter. 


Most  of  the  ISSAP  data  has  been  processed  for  sound  speed  at  frequencies  of  40,  65,  and  100  kHz. 

Our  initial  processing  efforts  have  concentrated  on  the  data  collected  in  Little  Bay  where  the  highest 
quality  backscatter  data  was  collected.  To  explore  possible  relationships  between  acoustic  backscatter 
and  sediment  properties,  the  EM3000  multibeam  backscatter  data  was  processed  with  a  technique 
called  Amplitude  Variation  with  Offset  (AVO)  [1]  a  technique  commonly  used  on  multi-channel 
seismic  reflection  data  for  the  exploration  and  characterization  of  subsurface  reservoirs.  AVO  analysis 
is  based  on  the  fundamental  assumption  that  the  seismic  amplitude  varies  with  the  offset  between  the 
seismic  source  and  detector  and  that  these  variations  are  inherently  related  to  the  acoustic  properties  of 
the  subsurface  reflectors.  Multibeam  sonars  acquire  backscatter  data  over  a  wide  range  of  incidence 
angles,  making  swath  data  well-suited  to  an  AVO  type  of  analysis. 

It  is  well  known  that  the  acoustic  backscatter  is  a  function  of  the  acoustic  properties  of  the  sediment, 
mainly  the  acoustic  impedance,  the  interface  roughness,  and  scattering  from  heterogeneities  within  the 
sediment  volume  0.  The  direct  determination  of  sediment  properties  from  remotely  measured 
backscatter  has  been  hampered  by  an  inability  to  separate  the  component  due  to  the  acoustic 
impedance  from  the  remainder  -  due  to  interface  roughness  and  volume  heterogeneities.  The 
application  of  AVO  analysis  to  multibeam  backscatter  attempts  to  address  this  ambiguity. 

Processing  of  the  multibeam  sonar  backscatter  data  described  here  is  applied  to  the  full-time  series 
data,  and  not  to  the  beam  averaged  backscatter  value.  The  AVO  analysis  is  applied  to  the  'best 
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estimate'  of  the  backscatter  due  to  seabed  reflection  -  that  is,  corrections  are  performed  to  remove  the 
time  varying  gain  (TVG)  applied  during  acquisition.  Additionally,  the  beam  angle  of  incidence  and 
ensonified  area  are  corrected  (during  data  acquisition  the  seafloor  is  assumed  to  be  flat)  based  on  the 
local  seafloor  slope,  a  correction  that  may  be  achieved  with  co-registered  multibeam  bathymetry  and 
backscatter. 

A  VO  is  applied  to  a  patch  of  seafloor  obtained  by  averaging  a  consecutive  number  of  pings,  or  angular 
response  (backscatter  as  a  function  of  angle  across  the  swath).  Processing  is  applied  independently  to 
the  port  and  starboard  angular  responses,  which  are  divided  into  three  range  intervals  based  on  angle. 
The  'near  range'  includes  grazing  angles  from  90°  to  65°,  the  'far  range'  from  65°  to  35°,  and  the  'outer 
range'  from  35°  to  5°.  Within  the  near  range,  the  mean  backscatter,  slope,  and  the  80°  intercept  of  the 
averaged  backscatter  angular  response  are  calculated  and  stored  as  AVO  attributes  (see  Fig.  3).  For 
the  far  range,  the  calculated  AVO  attributes  include  the  mean  backscatter,  slope,  and  the  intercept  at 
55°.  In  the  outer  range,  only  the  mean  backscatter  is  stored  as  an  attribute.  One  additional  AVO 
parameter,  the  'fluid  factor',  is  calculated  for  each  seafloor  patch.  For  that,  the  background  trend  line 
for  the  survey  is  defined  as  the  linear  regression  of  all  coordinate  pairs  (slope,  intercept)  in  the  slope- 
intercept  plane.  The  fluid  factor  attribute  is  defined  as  the  orthogonal  distance  of  each  coordinate  pair 
from  the  background  trend. 


Grazing  Angle  (Deg) 


Figure  3.  Backscatter  as  a  function  of  grazing  angle;  'A '  represents  the  near  range 
intercept  and  'B '  represents  the  near  range  slope. 


Once  the  AVO  parameters  are  calculated,  a  mathematical  backscatter  model  0  is  used  to  relate  the 
measured  backscatter  angular  response  of  each  seafloor  patch  to  seafloor  properties  of  acoustic 
impedance,  interface  roughness,  and  fluid  factor.  According  to  this  modeling,  the  fluid  factor  attribute 
is  directly  related  to  the  amount  of  free  fluid,  normally  gas,  in  the  sediment  structure.  The  same  AVO 
parameters  that  are  calculated  for  the  measured  backscatter  angular  response  are  calculated  for  the 
modeled  backscatter  angular  response.  A  model  inversion  is  performed  sequentially  by  adjusting  in 
order  the  near  range  slope;  near  range  intercept;  far  range  intercept,  far  range  slope,  and  the  fluid 
factor. 
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RESULTS 


Sound  speed  dispersion  was  observed  at  a  majority  of  stations  excluding  the  stations  with  a  sound 
speed  approximately  that  of  the  overlying  seawater.  Initial  comparisons  indicate  that  the  magnitude  of 
the  sound  speed  dispersion  is  related  to  the  mean  grain  size.  The  grain  size  analysis  conducted  by  URI 
was  initially  limited  to  five  grain  size  bins  which  did  not  provide  an  accurate  estimate  of  the  grain  size 
distribution.  URI  recently  completed  hydrometer  analysis  of  the  mud  component  (for  stations  with 
greater  than  10%  mud  fraction). 

The  AVO  approach  was  applied  to  the  Simrad  EM3000  multibeam  sonar  data  from  Little  Bay.  The 
estimated  impedance  and  was  compared  to  the  measured  ISSAP  sound  speed.  Shown  in  Figure  4  is  the 
estimated  impedance  ratio  as  a  function  of  the  ISSAP  measured  sound  speed  ratio  at  a  frequency  of  40 
kHz.  Initial  results  are  encouraging  and  indicate  a  strong  correlation  between  these  two  acoustic 
parameters.  Additional  field  work  is  required  to  include  a  larger  sample  of  sediment  types. 
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Figure  4.  Estimated  impedance  ratio  as  a  function  of  ISSAP 
sound  speed  ratio  at  40  kHz. 


IMPACT/APPLICATIONS 

The  ISSAP  has  provided  a  simple  and  quick  way  to  establish  the  lateral  distribution  of  sound  speed 
and  attenuation  variations  within  the  Geoclutter  area.  These  measurements  are  being  compared 
directly  to  backscatter  values  from  the  multibeam  system  and  to  the  predictions  of  impedance  and 
attenuation  made  from  the  Chirp  Sonar  by  Schock.  They  also  provide  information  on  the  range  of 
natural  variability  that  is  very  relevant  to  other  Navy  programs  (e.g.,  Capturing  Uncertainty  DRI  and 
Mine  Burial  Program  -  where  it  was  used). 

The  work  described  above  will  play  a  key  part  in  the  overall  development  of  robust  seafloor 
characterization  approaches,  particularly  through  helping  to  better  constrain  the  relationship  of  high- 
frequency  backscatter  to  seafloor  properties.  It  will  also  provide  critical  information  on  the 
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relationship  of  in  situ  properties  to  those  made  in  the  laboratory  as  well  as  those  extracted  remotely 
from  the  inversion  of  seismic  (chirp  sonar)  data. 

TRANSITIONS 

Data  requested  by  Chris  Jenkins  for  incorporation  into  global  sediment  property  database 

RELATED  PROJECTS 

Uncertainty  DRI,  Mine  Burial  DRI 
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